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I. Introduction 

Zero divided by zero is arguably the most important concept in calculus, as it is the 
gateway to the world of differentiation, as well as (via the fundamental theorem of 
^ calculus) the calculation of integrals. Organized methods of finding the right answers 

to zero-over-zero problems were developed by Newton and Leibniz and those who 
followed. The concept of limit was finally made rigorous long afterward by Cauchy 
and others. 

We are all sternly informed in grade school that zero over zero is impossible and 
must be shunned. A more accurate statement, however, (once we are ready to 
understand it) is that it is impossible to divide zero by zero out of context. 

All other arithmetic operations are independent of context, a rather remarkable 
fact. So six divided by three equals two, no matter whether you are sharing cookies 
with a couple of friends or making a fraction of a recipe of bread. Not so with zero 
divided by zero: if you take a photograph of two cars and observe that in the picture 
both cars travel zero feet in zero seconds, that is no help in comparing the cars' 
Cf) velocities. 

But if you can find some contextual information to a zero-over-zero problem and 
take a limit, you are back in business. As an example, if a problem traces back to 
the function 

| /(^<^ 

for x near 0, then we have at least four options for deciding what f(x) "ought to 
equal" when x = 0, namely plugging in nearby values and guessing, using algebra 
and cancellation, applying L'hospital's rule, or recognizing that the limit of / as x 
approaches zero represents the derivative of x 2 at 3. 

But what if the context of a zero-over-zero problem comes from a function of 
several variables? Is there a multivariable version of L'hospital's rule? As simple a 
limit as 

x-y 

lim 

(x,y) -> (0,0) sin a; - sin y 

x^y 

l 
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can be rather puzzling to consider, while something like 

7x 2 yz 5 + xy 3 — 3x 4 yz 

llHl 

(x,j/,z)-Ko,o,o) x 8 + x 2 y 2 z 4 + (y — x 3 + z 2 ) 2 + z e — xy 3 z 5 
is daunting indeed! 

A perusal of your multivariable calculus book will not likely help very much. At 
best you'll find some interesting individual examples worked out, with a caution not 
to try such things at home — at least not until you are properly warned against 
relying only upon individual lines through the origin! You can verify, for example, 
that if we set x = at, y = bt and z = ct, the above limit as t — > always works out 
to be zero, regardless of the values of a, b and c. But that is certainly not the whole 
story, as we shall demonstrate at the end of the paper. 

We seem largely to have overlooked the topic of finding an organized approach for 
resolving zero-over-zero limits of multivariable functions. 

The two papers |DSj and |Y], p. 71, both handle the specific situation of a 
two-variable indeterminate limit resolvable by taking the mixed second derivative 
d 2 jdxdy of the numerator and denominator functions. 

The paper |FKj has a version using first-order derivatives, but the theorem's use- 



fulness turns out to be limited, as we discuss after the proof of Theorem 2.1 below. 

There are some papers with a good treatment of the indeterminate limit of a 
quotient of a vector-valued function over a real-valued function, but these papers 
concern functions of a single variable. See [R], [A], [P] and [WJ for L'hopital -style 
theorems of this type. 

We now present a method for resolving many multivariable indeterminate limits. 

2. Nonisolated singular points 

A key aspect of a zero-over-zero limit of a multivariable function turns out to be 
the question of whether the singular point is isolated or not. Perhaps surprisingly, 
the nonisolated case turns out to yield the nicer theorem. 

If the singular point p is not an isolated zero of the denominator function g, then 
we understand the limit to be taken over points where g ^ 0. In order to hope for 
a limit to exist, the zero set of g must be contained in the zero set of / near p; 
otherwise the function will blow up near zeroes of g that are not zeroes of /. This 



hypothesis about the zero sets does not appear explicitly in Theorem |2.1[ but it is 
implied by the hypotheses that are there. 

Theorem 2.1 (Multivariable L'hopital's rule). Let M be a neighborhood in IR 2 con- 
taining a point p at which two differentiable functions f : J\f — > R and g : J\f — > R 
are zero. 
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Set 

C = {x G M : /(x) = g( X ) = 0}, 

and suppose that C is a smooth curve through p. 

Suppose there exists a vector v not tangent to C at p snc/i £/ia£ £/ie directional 
derivative D v g of g in the direction o/v is never zero within N ' . 

More generally, if C consists of a union of two or more smooth curves through p ; 
suppose that for each component Ei of Af\ C we can find a vector v i; not tangent at 
p to any of the curves comprising C , such that D v .g ^ on Ei. 

Then 

lim l^A= lim 
(*,v)-> P c/(x, y) ^ p A^fo y) 

i/ i/ie latter limits exist and are equal for all i. Each limit is assumed to be taken 
over the domain of points where the denominator is nonzero, and we assume in each 
case that p is a limit point of that domain. 

The same theorem holds for real-valued functions of n variables, with C a union 
of hyper surf aces. 

Proof. Suppose first that C is just a single smooth curve through p. Let L be the 
line through p in the direction of v. Since v is not tangent to C at p, parallel lines 
near L also cross through C near p; see Figure 1. 



L 




Figure 1. Sequences of points approaching a singularity 

Take a sequence (xi,yi) of points converging to p with g(xi,yi) ^ for all i. Let 
(a,, bi) be a corresponding sequence of points on C, also converging to p, such that 
(xi,yi) — (ai,bi) is a scalar multiple of v for all i. Now for each i, as in the proof 
of the single-variable L'hopital's rule, apply the Cauchy mean value theorem to the 
function f/g restricted to the line in the direction of v. Recalling that D v g ^ (it 
would be ironic to make a division-by-zero error in this proof!) we find that there is 
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a point (ci,di) in between (xi,yi) and (a,i,bi) such that 

Dvf(ci,dj) _ f(x l ,y i ) - fjaiJh) _ f(x h yj) 
D v g(ci,di) g{xi,yi) - gfabi) g(x h yi)' 

But now since the points (q, di) also converge to p, the desired result follows. 

Now if C consists of more than one curve through p then the same proof can be 
applied, with the domain restricted to any component Ei. The result is that the 
limit of f/g with domain Ei obeys the desired formula; since this is true (and gives 
the same value) for all of the finitely many components Ei, then the overall limit 
obeys the same formula. 

The proof for M. n goes through unchanged. □ 

We note here the importance of taking the directional derivative with respect to 
a fixed direction v. This works in the proof above because all lines passing near p 
in the direction of v do pass through the common zero set of / and g. A different 
strategy might be to take the directional derivative always in the direction straight 
out from a singular point. Indeed, the paper [FKJ, which handles only isolated 
singular points, does precisely this. The result is a true theorem, with proof similar 
to that above, but rarely is the expression obtained after differentiating any simpler 
than the original; the degree of a rational function, for example, does not decrease. 

Example 2.2. Let 

f(x,y) _ x-y 



h(x,y) 



g(x,y) sin x — sin y 

Then near the origin the zero sets of both numerator and denominator consist of the 
line y = x. Choose, say, v = (1,0), so that the directional derivative to be taken is 
just the partial derivative by x. Then 

lim = lim — = lim = 1. 

(x,y) -> (0,0) smx-smy g x cosx 

x^y 

Example 2.3. Let 

h(x,y,z) 



f(x,y,z) sinz — sin(x 2 + y 2 ^ 



g(x,y,z) tan(z — x 2 ) — tan(?/ 2 ) 

Then near the origin the zero sets of both numerator and denominator consist of 
the paraboloid z = x 2 + y 2 ■ Choose, v = (0, 0, 1), so that the directional derivative 



to be taken is just the partial derivative by z. The hypotheses of Theorem 2.1 are 
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satisfied, and 



sin z — sin(x 2 + y 2 ) cos z 

lim ; — - = hm 



(x,i/,z)-Ko,o,o) tan(z — x 2 ) — tan(y 2 ) (x,j/,z)-»-(o,o,o) sec 2 (2 — a; 2 ) 
Example 2.4. Let 

g(x,y) cos x — cosy 
Near the origin, numerator and denominator are zero when y = ±x. Choosing 
V! = (1,0), we note that D VL g = — sinx is zero on the y axis, which intersects 
neither the open "east component" E\ of the plane bounded by the northeast and 
southeast arms of the zero set of g, nor the opposite "west component." We calculate 

2x 

lim h(x, y) = lim = —2 

(s,iO-Ko,o) (x,y)->(o,o) sinx 

for (x,y) restricted to E\. 

For the north and the south components we use v 2 = (0, 1) and obtain the same 
limit, so the overall limit is —2. 

3. Isolated singular points 

Now we move to the case where the denominator function g is (or approaches) zero 
at an isolated point p. We prescribe the following algorithm, not always conclusive, 
for deciding the value of such a limit. The basic idea is to use a change of variables 
that makes the denominator into a sum of squares. 

Algorithm: To attempt to resolve a limit in W 1 of the form 

,. /(*) 
hm — — 

where g(x) ^ on a deleted neighborhood of p but g(x) equals (or approaches) at 
the point p: 

Step 1: Make a simple, preliminary substitution, such as holding all but one variable 
constant, and find the resulting restricted limit at p. 

• If the restricted limit equals 0, proceed to Step 2. 

• If the restricted limit does not exist, then of course the unrestricted limit 
does not exist, either. 

• If the restricted limit exists and equals I 7^ 0, replace the numerator /(x) 
by /(x) — £g(x), so that the restricted limit will now equal zero. 

Thus, the two remaining possibilities are that the (unrestricted) limit is 
zero or undefined. 
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Step 2: See whether a separation technique such as in Example 3^ below resolves the 
limit. If not, proceed to Step 3. 

Step 3: Write the denominator g(x) as a sum of squares «i(x) 2 + • • • + m„(x) 2 plus 
(optionally) another function v (x). Each itj(x) must approach zero as x — > p. 

The function v, if present, will need to be less significant near p than the 
sum of Mj(x) 2 , as described further below. 

Step 4: To prove the limit does not exist, try the substitution Mj(x) = mit for all i. If 
this defines, for each appropriate choice of constants {m^}, a curve through 
p, and if the resulting restricted limit does not exist or depends on the values 
of {rrii}, then the original limit does not exist. 

Step 5: To prove the limit is zero, make a polar or spherical coordinate substitution 
iij(x) = hi(8i, . . . , O n -i)p, where the sum of squares of the functions hi is 
identically 1. If the numerator /(x) and the residual denominator term v(x) 
can both be written as linear combinations of p aj times bounded functions, 
with the powers <x,- all strictly greater than 2, then we can divide top and 
bottom by p 2 , leaving a positive power of p times a bounded function. Since 
p — > as x — > 0, the overall limit is zero. 

Example 3.1. We begin with an example that can be resolved by separating vari- 
ables: 

2 i -4 

,. x + sm y 
hm — 5 -. 

(x,i/)-Ko,Q) sin x + y 4 

Setting, for example, y — and applying the single-variable L'hospital's rule we get 
a restricted limit of 1. Thus, according to Step 1 we should try to prove that the 
following new limit equals zero: 



lim 



x 2 + sin 4 y — sin 2 x — y A 



Now this equals 



(s,iO->-(o,o) sin x + y q 



x 2 — sin 2 x sin 4 y — y 4, 

hm — t, 7 + hm 



(cc,y)->(o,o) sm x + y* O,y)-Ko,o) sm x + y q 

which is bounded in absolute value by 

\x 2 — sin 2 a;| I — sin 4 y + y 4 \ 

hm ^ h hm 



(x,y)^(o,o) sin 2 x (x,y)-+(o,Q) y 4 
L'hospital's rule shows each of these limits to be zero, as desired. 

Example 3.2. We continue with a simple example, standard in textbooks: 

f(x,y) _ xy 
g(x,y) x 2 + y 2 ' 
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The singularity at the origin is isolated, and the preliminary substitution y = gives 
a restricted limit of 0. 

The denominator is already a sum of squares of polynomials, with U\{x,y) = x 
and u 2 (x,y) = y. As prescribed in step 3 of the algorithm, we set 

x = m\t and y = m 2 t, 

resulting in the limit 

m\m2t 2 17111712 
lim = ' 

t->o m\t 2 + m\t 2 m\ + m\ ' 

since this limit depends on m\ and m2, the original limit does not exist. 

The point commonly made in connection with the above function is that although 
you can prove that a limit does not exist by restricting to straight lines through the 
origin in the domain, this method does not suffice to prove that a limit does exist. 
A typical example cited is the following, which we will analyze first by lines through 
the origin and then through the algorithm above. 

Example 3.3. Consider 

v x 2 y 

hm — j -. 

(x,y)-*(0,0) x 4 + y 2 
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when r / 0. Holding 9 constant and taking the limit as r — > we always get a 
(restricted) limit of 0. This is reflected by the lighter curves in the figure (on which 
9 is constant) all converging to one center point. 

We might think, then, that the overall limit of g is zero. As the textbooks point 
out, this is not the shown by the substitution y = x 2 . Our algorithm 

prescribes essentially the same substitution; with 

u 1 (x,y)=x 2 and u 2 (x,y)=y, 

we set 

x 2 = mit and y = m 2 t. 

The limit as t — )■ of the restricted function is now 

mim 2 
m\ + m\ 

This is different for different values of m, so the original limit does not exist. This is 
reflected by the heavier curves in the figure (on which r is constant) having to drape 
up over the hump. 



Example 3.4. Let 



x 3 y 3 

h(x,y) 



x 6 + y A ' 
We set 

x 3 = mit and y 2 = m 2 t. 

Then 

_ m^ff 5 / 2 m x m 3 2 12 t i/2 
(ml + ml)t 2 (ml + m|) 

which goes to zero at t — > 0. Suspecting that the limit is zero, we move to step 4. 
Set x 3 = r cos 9 and y 2 = r sin 9. Then 

r 5/2 | cos 6*1 1 sin^l 3 / 2 1/9 
\h\ = ^ !— < r 1 / 2 , 



so that the limit is zero. 



Note: The substitution above does cover all values of x and y, whether positive 
or negative. The substitution y 2 = r sin 9 is not one-to-one; it has to use the same 
values of 9 & [0, ir] for positive and negative values of y. But since the polar limit is 
zero and since \\(x, y)\\ small implies r is small, the original limit is also zero. 



A fancier example features more than two square terms in the denominator. 
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Example 3.5. Let 



h(x,y) 



S 2 

x y 



f{x,y) _ 

g(x,y) x 6 + x 2 y 2 + y e 




Figure 3. z — x 3 y /(x b + x y + y b ) 

The algorithm asks for just two square terms from the denominator. We try 
different pairs. 
First let us set 

u 1 (x,y)=x 3 , u 2 (x,y) = y 3 , and v(x,y) = x 2 y 2 . 

Then (setting the mi's equal to 1 for quicker analysis) we make the substitution 



and obtain 



x 3 = t = y 3 



t 5/3 



t l/3 



t 2 + t 4 / 3 + t 2 l + 2t 2 /3' 

which goes to zero. But if we try to verify the zero limit by setting x 3 = r cos 9 and 
y 3 = r sin 9 we get 

r 5 cos 3 9 sin 2 9 
r 2 + r 4 / 3 (cos#sinfl) 2 / 3 ' 

which is inconclusive. Note that with this substitution, the power of t coming from 
the side term v(x,y) was smaller than 2. 
Next we try 

ui(x,y) = x 3 , u 2 (x,y) = xy, and v(x,y) = y 6 . 
The substitution 

x 3 = mit and xy = m 2 t 
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gives 



{x)(xy) 2 (x)(xyY 



x & + x 2 y 2 + y 6 x 6 + x 2 y 2 
m\t 2 + m\t 2 ' 

which is bounded by a constant times a positive power of t, and so goes to zero. We 
verify this zero limit by setting 

a: 3 = rcos#, xy = r sin 9 
, . r 7 / 3 | cos 1 / 3 6 1 sin 2 #1 r 7 / 3 



< = r 



1/3 



^•2 j 

The limit of r 1//3 is zero, so that the original function has a zero limit at the origin. 

One final example illustrates how we may replace functions by Taylor polynomial 
approximations, and also use dominant square terms to eliminate non-dominant 
terms. 

Example 3.6. Let 

2-2cos(xV) 



h(x,y) 



x 10 + x e y 2 + y 6 — x 9 sin y 
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Define 



9i(x,y) 


= x 10 + x 6 y 2 + y 6 


92(x,y) 


= —x 9 smy 


9~2(x,y) 


= -x y 


f(x,y) 


= 2 - 2 cos (a V) 


f(x,y) 


4 4 

= x y . 



Now by setting u = x y and using the single-variable L'hopital's rule we see that 



similarly we verify 



lim |H = l; 



lim ^4 = 1. 

(x.yj-Ko.o) g 2 {x,y) 



Then if we show that the limits of §ijg\ and f/g 1 are zero we will also know that 
the limit of gijg\ is zero, and the limit of (gi + g-zj/gi is 1; thus 

r u r f r f 9i + 92 f ,. / n 
lim h = lim = lim = lim — = 0. 

gi + g2 gi + g2 gi i gi 

To analyze g~2/gi, set 

x 3 y = r cos 9 and x 5 = r sin 9. 

Then 

9~2 {x, ?/) r| cos^| r 6//5 | sin^l 6 / 5 
Si(a;,y) r 2 + y< 

< 

which goes to zero as r — > 0. Thus, 



6 



11/5 

< — = 



lim ^4 = . 

(*,y)-Ko,o) gi{x,y) 

To analyze //<?i, this time we need to highlight a different pair of square terms 
from the denominator; set 

x 3 y = r cos 9 and y 3 = r sin 9. 

Then 

I f(x,y) | _ r 4 / 3 |cos#| 4 / 3 r 8 / 9 |sin#| 8 / 9 
y) r 2 + x 10 
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,.20/9 

< — = r ' /9 ^ °- 

One technicality is that some of the limits being combined here, such as the limit 
of ///, are undefined when x = or y = 0. For this, note that if the overall limit 
of h is zero when restricted to each of finitely many domains then the limit over the 
union of those domains is also zero. Clearly the limit of h is zero when restricted to 
x = or to y = 0, so the analysis is complete. 

4. Unfinished business 

We promised in the introduction to complete the analysis of the limit 

7x 2 yz 5 + xy 3 — 3x A yz 

hm ■ — -. 

(x,y,z)-*(o,o,Q) x 8 + x 2 y 2 z 4 + (y — x 3 + z 2 ) 2 + z 6 — xy 3 z 5 

We noted there that all restricted limits along lines through the origin are zero. To 
analyze further, set 

4 

U\ = X 

u 2 = xyz 2 

u 3 = y — x 3 + z 2 

m 4 = z 3 

we obtain 

lim ^ m i /4m2M4 + u i /4 ( u 3 + v 3 / 4 - u 2 / 3 ) 3 - 3ui(m 3 + u^ 4 - u/ 3 )uj 3 

u\ + u\ + u\ + u\ - u 2 u A (u 3 + u 3 / 4 - uj 3 ) 2 

Now the non-square denominator term and all numerator terms but one are of degree 
greater than 2. But the last numerator term is of degree 2. 

The substitution in Step 4 of the algorithm of section [3] would result in an overde- 
termined system. Since the term of degree 2 in the numerator only involves u±, u% 
and «4, we leave out u 2 in determining t. Also, to eliminate fractional exponents we 
will use t 12 in place of t. Choosing mi = = = 1 we set t 12 = u\ = U3 = M4 
and solve for x, y, and z to obtain 

x = t 3 

z = t 4 

y = t 12 + t 9 - t S . 

Then the original limit translates to 

3t 24 + 0(t 25 ) 
3t 24 + (t 25 ) ' 
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whose limit is 1. So the original function is bounded in a neighborhood of (0,0,0) 
(since no numerator term has degree less than 2), and it turns out that the function 
goes to zero except on a tiny thread through the origin. 
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